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Abstract 

We report on a search for heavy neutrinos in .B-meson decays. The results are obtained using 
a data sample that contains 772 x 10 6 BB pairs collected at the T(45) resonance with the Belle 
detector at the KEKB asymmetric-energy e + e~ collider. No signal is observed and upper limits 
are set on mixing of heavy neutrinos with left-handed neutrinos of the Standard Model in the mass 
range 0.5 GeV/c 2 - 5.0 GeV/c 2 . 

PACS numbers: 12.60.-i,13.35.Hb,14.60.Pq 
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The masses of particles in the Standard Model (SM) are generated by the coupling of the 
Higgs field to the left- and right-handed components of a given particle. There being no right- 
handed neutrino components in the SM, neutrinos in the SM are strictly massless. However, 
experimental data on neutrino oscillations show that neutrinos are not massless, though their 
masses are very small [l[. Therefore, a mechanism beyond the SM is needed to establish 
neutrino masses. One possibility is the addition of right-handed neutrinos, which may also 
have a Majorana mass, naturally explaining the smallness of the observed neutrino masses 
via the so-called "see-saw" mechanism [2|. For example, the neutrino minimal Standard 
Model (VMSM) [3] introduces three right-handed singlet heavy neutrinos, so that every left- 
handed particle has a right-handed counterpart. This model explains neutrino oscillations, 
the existence of dark matter and baryogenesis with the same set of parameters. Heavy 
neutrinos also appear in other extensions to the SM, such as SUSY [4|, grand unification 
theories [E| or models with exotic Higgs representations [6[. 

In general, neutrino flavor eigenstates need not coincide with the mass eigenstates but 
may be related through a unitary transformation, similar to the one that applies to the 
quark sector, 

Vg. = ^ UgjVi, ol = e,/z,T, 1 = 1,2,3,4,... (1) 

i 

where a denotes the flavor eigenstates and i denotes the mass eigenstates. Production 
and decay diagrams for heavy neutrinos are shown in Fig. [IJ The coupling of the heavy 
neutrino v± to the charged current of flavor a is characterized by a quantity U a ±. Below, 
we denote a heavy neutrino in the mass range accessible at Belle and its corresponding 
coupling constant by Vh and U a , respectively. Existing experimental results are reviewed 
and discussed in Ref. p]]. 

In this paper, we describe a direct search for heavy neutrino decays v h — » £ ± 7r =F , I — e, // 
with the Belle detector. The measurement is based on a data sample that contains 772 
million BB pairs, which corresponds to 711 fb -1 , collected at the T(AS) resonance with 
the Belle detector operating at the KEKB asymmetric-energy e + e~ collider [8j. The Belle 
detector is a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold Cherenkov 
counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located inside a 
superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux return 
located outside the coil is instrumented to detect K\ mesons and to identify muons (KLM). 
The detector is described in detail elsewhere [9|. Tracking at Belle is done using the SVD and 
CDC. Two configurations of the SVD, referred to as SVD1 and SVD2, were used; 140 fbT 1 
were accumulated with SVD1 and the remainder with SVD2. SVD1 consists of three layers 
of double-sided silicon strip detectors arranged in a barrel surrounding a 20 mm beam pipe 
at radii of 30 mm, 45.5 mm and 60.5 mm. SVD2 has a 4-layer configuration with radii of 
20 mm, 43.5 mm, 70 mm and 80 mm surrounding a 15 mm radius beam pipe. The inner and 
outer radii of the CDC are 83 mm and 880 mm. 

Backgrounds are studied using Monte Carlo (MC) samples of known BB decays from 
b — > c processes (generic MC) that have three times the statistics of the Belle dataset. Sig- 
nal MC samples of 500,000 events each for different heavy neutrino masses and production 
mechanisms are used to evaluate the response of the detector, determine its acceptance 
and efficiency, and optimize selection criteria. Events are generated using the EvtGen pro- 



gram 10]. Heavy neutrinos are produced and decayed using a phase space model. 
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At Belle, the most favorable mass range to look for a heavy neutrino is M(K) < Miv^) < 
M(B) [ll|. This analysis uses the leptonic and semileptonic B meson decays B — > X£uh, 
where l — e y \i and X may be a charm meson £)W, a light meson (jr, p, i], etc.) or 'nothing' 
(purely leptonic decay), with relative rates as given in Ref. 

A distinctive feature of the heavy neutrino is its long expected flight length: for Miv^) = 
lGeV/c 2 and \U e \ 2 = |£7 M | 2 = 10~ 4 the flight length is cr ~ 20m. Therefore, the expected 
overall reconstruction efficiency is small. To improve sensitivity, a partial reconstruction 
technique is used. A candidate is formed from two leptons and a pion (£2^1^), where £\ and 
7r have opposite charge and form the heavy neutrino candidate with a vertex displaced from 
the interaction point (IP). The lepton l\ is referred to as the 'signal lepton,' while the lepton 
£2, which comes from the B decay, is referred to as the 'production lepton.' In this analysis, 
the heavy neutrino is assumed to be a Majorana fermion and may decay to a lepton of any 
charge regardless of the original 5-meson flavor. If the heavy neutrino were a Dirac fermion, 
the production and decay leptons would necessarily have opposite charge. 

If the heavy neutrino is light enough to be produced via B — > D^Hu^, these production 
modes are expected to dominate over decays to light mesons due to the small value of the 
ratio of the relevant CKM matrix elements |14&|/|V^|. The background is more severe for 
smaller heavy neutrino masses, M(z/^) < 2GeV/c 2 , so an analysis using only B — > D^£vh 
modes is used this "small mass" regime, while the full inclusive analysis is used in the "large 
mass" regime. 

To suppress the QED background, the charged multiplicity in the event is required to be 
larger than four. Charged tracks positively identified as electrons or muons (as defined in 
the next paragraph) with laboratory-frame momentum greater than 0.5GeV/c are used as 
leptons. All other tracks in the event are treated as pion candidates. Additional selection 
criteria for the lepton and pion tracks are described below. 

A significant background remains for heavy neutrino candidates from particles with sim- 
ilar event topology, notably Kg — > 7r + 7r~, A — > pn~ , 7 — > e + e~. These backgrounds are 
suppressed by strict lepton identification requirements. Electrons are identified using the 
energy and shower profile in the ECL, the light yield in the ACC and the specific ionization 
energy loss in the CDC (dE/dx). This information is used to form an electron (£ e ) and non- 
electron (£ g ) likelihood; these are combined into a likelihood ratio TZ e = C e /(C e + C s ) fl2j ]. 
Applying a requirement on lZ e , electrons are selected with an efficiency and a misidentifica- 
tion rate of approximately 90% and 0.1%, respectively, in the kinematic region of interest. 
Muons are distinguished from other charged tracks by their range and hit profiles in the 
KLM. This information is utilized in a likelihood ratio approach [13j similar to the one used 
for the electron identification. Muons are selected with an efficiency and a misidentification 
rate of approximately 90% and 1%, respectively, in the kinematic region of interest. These 
requirements are reversed in order to produce a lepton veto for identifying pion candidates. 

We select well-vertexed heavy neutrino candidates using dr, the distance of closest ap- 
proach to the IP in the plane perpendicular to the beam axis for each track; d<p, the angle 
between the momentum vector and decay vertex vector of the heavy neutrino candidate; and 
dz vtx , the distance between the daughter tracks at their closest approach in the direction 
parallel to the beam. Requirements vary depending on the presence of SVD hits on the 
tracks and on the heavy neutrino candidate flight length. For vertices closest to the IP that 
include tracks with SVD hits, we require dr > 0.09 cm, dcp < 0.03 rad and dz vtx < 0.4 cm, 
while for farthest vertices, more farther than 30 cm from the IP that do not include SVD 
hits, the requirements are dr > 5 cm, d(p < 0.15 rad and dz vtx < 2 cm. The efficiency of 
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these requirements for background is 0.8%, while for signal it is 50-80%, depending on the 
mass and production mode of the heavy neutrino. The signal lepton and pion are fit to a 
common vertex. Only candidates with x\/ndf < 16, where Xi is the goodness of fit and ndf 
is the number of degrees of freedom, are accepted. A second vertex fit of the heavy neutrino 
candidate and the production lepton is performed with the vertex constrained to the IP; 
candidates with \\jndj < 4 are retained. 

For combinatorial background, the daughter tracks of the heavy neutrino candidate often 
originate from the vicinity of the IP rather than the candidate's decay vertex. In order 
to suppress this background, the difference between the radial coordinates of the closest 
associated hit in the SVD or CDC of either of the two daughter tracks to the IP (rg or r n ) 
and the candidate's decay vertex (r vtx ) is calculated: dr^ = min(r^,r 7r ) — r vtx . The analysis 
requires dr^ > —2 cm. 

For the "small mass" (M(u h ) < 2 GeV/c 2 ) analysis, B -> D&lv h events are selected 
using the recoil mass against the ££ir system. This requirement is related to the kinematics 
of the decay under study. For B — ► X£uh — > X££n decays, the mass of X can be obtained 
from M| = (EcM—Eun) 2 — Pe£ n —PB~^Pu7T - PB, where Eqm and Pb are the B meson center- 
of-mass (CM) energy and momentum and Eg^ and P^ are the CM energy and momentum 
of the ££tt system. The last term in this equation cannot be calculated as the B direction 
remains unknown, so we redefine the recoil mass as M\ = (Ecu ~ Eu n ) 2 — P 2 £w — P|. For 
events with X = D^*\ the Mx distribution has overlapping peaks around the masses of the 
D and D*, while for background events the recoil mass has a broader distribution, as shown 
in Fig. [2J Events with 1.4GeV/c 2 < M x < 2.4GeV/c 2 are selected as candidates. 

To reject protons from the decays of long-lived baryons, we impose a loose proton veto 
for the pion candidate. For each track, the likelihood values C p and Ck of the proton 
and kaon hypotheses, respectively, are determined from the information provided by the 
hadron identification system (CDC, ACC, and TOF). A track is identified as a proton if 
Cp/(C P + Ck) > 0.99. Background events, rejected by the veto, are concentrated at heavy 
neutrino masses below 2GeV/c 2 and thus this veto is applied in the "small mass" analysis 
only. 

Using the requirements described above, the number of background events is reduced by 
a factor of ~ 10 6 to a handful of events, as shown in Fig. [3j 

If the heavy neutrino lifetime is long enough, then the number of neutrinos detected in 
the Belle detector is (in units where h = c = 1) 



where Nbb is the number of BB pairs, BiB — > X£v^) is the total branching fraction for 
u h production, B(u h — > £n) is the branching fraction of the reconstructed decay, e(R) is the 
reconstruction efficiency of the Uh decaying at a distance R from the IP and m, p and T are 
the mass, momentum and full width of the heavy neutrino, respectively. Additionally, to 
factor out the \U\ 2 dependence, we define \U a \ 2 fi(m) = B[B — > X£v^) and \U ^ j2{jn) = 
r(z^i — > in) = B{i>h — > £n)T, where a and (3 denote the flavor of the charged lepton 
produced in the B and v h decay, respectively. The exponent in the integrand of Eqn. (T5]) is 
approximated by unity. Integration is performed over the full volume used to reconstruct the 




dR 



(2) 



6 



heavy neutrino vertex, which depends on the reconstruction requirements. The expressions 
for BiB — > Xlvh) and T(vh — > £tt) are taken from Ref. [1 1| , and require only very general 
assumptions (i.e., they are not specific to z/MSM). 

The reconstruction efficiency e(R) for signal events does not depend significantly on the 
reconstruction mode (ee7r, /i/i7r or e/i7r), but does depend strongly on the heavy neutrino 
mass. The maximum efficiency is achieved at M(z//J ~ 2GeV/c 2 . For a given mass, the 
efficiency also depends on the 5-meson decay mode in which the heavy neutrino is produced. 
The efficiency for leptonic decays is higher than for semileptonic decays, and the efficiency 
for B — y iriuh decays is higher than for B — > Dlv^ decays. For example, for a heavy neutrino 
produced in B — > Dtvh with M{i>h) = 0.5 GeV/c 2 (2 GeV/c 2 ), the average reconstruction ef- 
ficiency is 3.1% (3.8%); for those produced in B ir£v h with M(y h ) = 2 GeV/c 2 (4 GeV/c 2 ), 
the efficiency is 7.7% (1.3%); and for those with M(z//J = 3.0 GeV/c 2 produced in B — > tv h 
decay, the efficiency is 10%. The efficiency drops with the radius R of the decay vertex from 
the beam axis. The effective range of neutrino reconstruction extends to R ~ 60 cm. 

The calculated total branching fraction for heavy neutrino production B(B — > Xtvy/) for 
the "small mass" and "large mass" analyses corresponds to 

B(B -> Xlv h ) sman mass = B(B -> D£u h ) + B(B -> D*£u h ) (3) 

and 

B[B -> Xiv h ) 

large mass 

i 

respectively, where the summation is done over D, D*, 7T, p, rj, rj', u, <f> and 'nothing.' These 
are not exact expressions but rather estimates of lower bounds on BiB — > Xivh), which lead 
to conservative upper limits on \U\ 2 . 

The systematic uncertainty of each of the event selection criteria is estimated from the 
difference in the efficiencies obtained in data and MC. A summary of all systematic uncer- 
tainties is presented in Table 0] Since all particles used in the systematic uncertainty study 
decay relatively close to the IP compared to the expectation for a heavy neutrino, we require 
where possible that the decay vertices be farther than 4 cm from the IP in the transverse 
plane to put more weight on large decay lengths. To estimate the systematic uncertainty 
due to tracking, we compare the number of fully and partially reconstructed D* decays in 
the decay chain D* — > Dtt + , D — > Kgirir, Kg — > tttt, where in the latter case one of the 
pions from the Kg is explicitly left unreconstructed. To estimate the systematic uncertainty 
of the recoil mass requirement, we reconstruct B — > DDi*\ D — > K s ir events and study 
the mass recoiling against the D-meson. The D decay topology is similar to here, and 
we treat the difference in recoil mass efficiency between data and MC as the systematic 
uncertainty of the recoil mass requirement. To estimate the systematic uncertainty of the 
electron identification, we reconstruct tt° — > 77 events, where one of the photons converts 
into e + e~ in the detector and one of these conversion particles is identified as an electron. 
The difference of the identification efficiency of the other daughter between data and MC 
is treated as a systematic uncertainty. For the muon identification, we perform a similar 
study with a J /if} — > fi + fi~ sample. To estimate the systematic uncertainty of other recon- 
struction requirements, we apply these requirements to K s decays, which have a topology 
similar to heavy neutrino decays. All systematic uncertainties are summed in quadrature, 
leading to total systematic uncertainties of 31.9% and 31.6% for the "small mass" and "large 
mass" regimes, respectively. The largest contributions to the systematic uncertainties are 
X2 (10.1%), dz vtx (10.0%) and tracking (8.7% per track, added linearly). 
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TABLE I. Summary of systematic uncertainties. Asterisks indicate requirements used in the "small 
mass" analysis only. 
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After all the event selection criteria were fixed from the MC study, the data were analyzed 
and the coupling constants \U e \ 2 , |t/ M | 2 and |LV e ||f/ M | were obtained separately using the decay 
modes eeii, fifiir and e/i7r + [ten, respectively. Distributions of the heavy neutrino mass in 
generic MC and data are shown in Fig. [31 In agreement with MC expectations, only a few 
isolated events are observed and we set upper limits on \U\ 2 according to Ref. taking 
into account the systematic uncertainty calculated above. For non-empty bins and empty 
bins far from non-empty bins, we set Poisson upper limits, assuming small background, as 
suggested from the MC study. In the vicinity of non-empty bins, we use Gaussian fits to 
interpolate between empty and non-empty regions. The widths of the Gaussians are fixed 
from MC. We use bins of 3 MeV/c 2 width, since the mass resolution evolves from ~ 3 MeV/c 2 
at M(v h ) = 1 GeV/c 2 to ~ 12MeV/c 2 at M(u h ) = 4GeV/c 2 . The resulting upper limits at 
90% CL on the number of events and coupling constants are shown in Fig. HJ 

In conclusion, upper limits on the mixing of heavy right-handed neutrinos with the con- 
ventional SM left-handed neutrinos in the mass range 0.5 — 5.0 GeV/c 2 have been obtained. 
The maximum sensitivities are achieved around 2 GeV/c 2 and are 3.1 x 10~ 5 , 3.1 x 10~ 5 and 
2.2 x 10 -5 for \U e \ 2 , |t/ M | 2 and |C/ e ||C/ M |, respectively. The corresponding upper limit for the 
product branching fraction is B(B — > Xlvh) x B{i>h — > £n + ) < 7.5 x 10 -7 for £ = e or ji. A 
comparison with existing results for \U e \ 2 and \U^\ 2 is shown in Fig. [5l 
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FIG. 3. Distributions of M{v^) for ee7r, /U^7r and efin + fieir reconstruction modes in generic MC 
(unsealed) (a), and data (b). The dotted line shows the boundary between the "small mass" 
and "large mass" methods. The filled (black) histograms are for candidates with opposite-charge 
leptons, while the open (white) histograms are for candidates with same-charge leptons. 
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FIG. 4. Upper limits at 90% CL on the number of signal events (a) and \U e \ 2 , \U^\ 2 and 1^11^1 
(b). The dotted line shows the boundary between the "small mass" and "large mass" methods. 
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